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A quasi- one- dimensional quantum dot containing two in- 
teracting electrons is analyzed in search of signatures 
of chaos. The two- electron energy spectrum is obtained 
by diagonalization of the Hamiltonian including the ex- 
act Coulomb interaction. We find that the level-spacing 
fluctuations follow closely a Wigner-Dyson distribution, 
which indicates the emergence of quantum signatures of 
chaos due to the Coulomb interaction in an otherwise 
non-chaotic system. In general, the Poincare maps of 
a classical analog of this quantum mechanical problem 
can exhibit a mixed classical dynamics. However, for the 
range of energies involved in the present system, the dy- 
namics is strongly chaotic, aside from small regular re- 
gions. The system we study models a realistic semicon- 
ductor nanostructure, with electronic parameters typical 
of gallium arsenide. 

PACS: 73.23.-b, 73.61.-r, 05.45. G 



In a many-body system, it is possible for signatures 
of quantum chaos to appear due solely to the interac- 
tions among its particles. During the last decade, such 
interaction-induced signatures of quantum chaos havj 
been investigated, for example, in spin-fermion models,! 
in the compound nuclear state with 12 particles in thft 
sd-shellH and in the heavy rare-earth atom of Cerium.El 
In those studies, the evidence for quantum signatures of 
chaos found in the level-spacing statistics or in the statis- 
tical properties of the wave functions was considered to 
be conclusive. This conclusion is not entirely surprising 
since all of those systems, having relatively large numbers 
of particles (> 10), are similar to the complex nuclear 
systems that motivated the introduction of the-|ideas of 
random matrix theory (RMT) in the first placeD On the 
other hand it is not obvious a priori, whether the inclu- 
sion of the interaction in few-body systems, like for ex- 
ample currently available semiconductor quantum-dots, 
leads also to signatures of quantum chaos. In a recent 
study, the effect of electronic interactions was consick 
ered in a parabolically-confined three-electron system.El 
(It is well known that three is the lowest number of in- 
teracting electrons necessary to break the integrability in 
a parabolic quantum dot.) In that study, the crossover 
from regular to irregular spectra as a function of the in- 
teraction strength was found to be incomplete, possibly 
due to the existence of hidden symmetries not taken into 
account in the statistical analysis. Therefore, the ques- 
tion of the emergence of signatures of quantum chaos due 



to interparticle interactions in systems of very few parti- 
cles remains open. 

Closely related to the issue of characterizing the dy- 
namical properties of simple interacting systems is the 
problem of quantum control with external fields. The 
manipulation of few particles (electrons and holes) in 
semiconductor quantum dots is a potentially impor- 
tant techuological problem that is receiving increasing 
attentionJj'El In this context, recent theoretical studies 
have shown interesting effects of single-electron turnstile 
behavior, and localization and correlations in systenia of 
quantum dots with two interacting electrons in them.EI In 
the present Letter we investigate the signatures of quan- 
tum chaos in a similar system, ie., two interacting elec- 
trons in a quasi-one-dimensional semiconductor quantum 
dot. We show that the Coulomb interaction between 
the electrons induces an unambiguos transition from a 
regular spectrum to a spectrum that follows closely the 
predictions of RMT for systems whose classical analog 
exhibits chaotic dynamics. 

In order to fully characterize the emergence of chaos 
due to interactions in this simple system, we also study 
the dynamics of its classical analog. The Poincare maps 
that describe the classical system show a strongly (albeit 
mixed) chaotic behavior due to the inclusion of the true 
Coulomb interaction in the system. 

We assume that the quantum dot has narrow parabolic 
confinement in the transversal x — y dimensions, so that 
the energies associated with those modes are high com- 
pared to the energies of the remaining degree of freedom 
(Born approximation). The two-electron wave function 
can then be written as 



*(ri, ra) = (t>(xi)(l){yi)(l){x2)(t>{y2)'^{zi, Z2), 



(1) 



where (j){x) is the lowest harmonic oscillator energy eigen- 
state. The energy eigenstates satisfy 

[Hq{zi) -f Ha{z2) + Vid{\zi - Z2|)]$(zi, Z2) 

^E^{z^,Z2), (2) 

where Ho{z) = -{h^ /2m*) d'^jdz^ + V{z) is the single- 
particle Hamiltonian with V"(z) being the quantum-dot 
defining potential, ra* is the effective mass, and Y\d is 
the Coulomb interaction given by 



y\D[\z\ - Z2\ 



dx\dyidx2dy2 



e^^\xi)cbHyiW{^2W{y2 

e|ri - r2| 



(3) 



We use the values of the dielectric constant e — 13.1 and 
m* = 0.067?ne corresponding to gallium arsenide. We 
choose to work with a quasi-one-dimensional semicon- 
ductor quantum dot confined in 15 A in the x-y plane 
and a width of 800 A in the z-direction. 

In the absence of the interaction term in Eq. (||), 
the Hamiltonian is a sum of two single-particle one- 
dimensional Hamiltonians, whose classical counterpart 



is obviously integrable. The main question we seek to 
answer is whether the Coulomb interaction between the 
electrons introduces chaos in the system. 

In order to look for signatures of quantum chaos, we 
follow a standard statistical analysis of the energy spec- 
trum, which consists of the following steps. First we cal- 
culate the exact spectrum {En} by diagonalization of the 
Hamiltonian matrix. The level spectrum is used to ob- 
tain the smoothed counting function Nav{E) which gives 
the cumulative number of states below an energy E. In 
order to analize the structure of the level fluctuations 
properties one "unfolds" the spectrum by applying the 
well kwown transformation Xn = Nav{En).u From the 
unfolded spectrum one calculates the nearest-neighbor 
spacing (NNS) distribution P{s), where Si = Xi^i — Xi is 
the NNS. 

We first consider the spectral properties of the non- 
interacting two-electron problem whose Hamiltonian is 
Ho{zi) + Ho{z2)- Its eigenstates can be classified by their 
total spin in singlets and triplets. To compute the NNS 
distribution we use eigenstates of a given spin. In the in- 
set of Fig. 1 we show the obtained NNS non-interacting 
distribution P/v/(s) (histogram) which follows a Poisson 
distribution (characteristic of an uncorrelated sequence 
of energy levels) given by Pp{s) ~ e"" and shown for 
comparison as a solid thin line. Due to the finite di- 
mension of the Hilbcrt space, Nav{E) saturates in the 
highest energy region. Therefore, we compute the-^NNS 
distribution using the lowest ~ 1000 eigenvalues. t3 The 
obtained Poisson distribution is an expected signature of 
most quantum two-dimensipaal systems whose classical 
counterparts are integrableJiJ 

To analize the interacting spectrum we diagonalize ex- 
actly Eq. (|2|). We also take into account the symmetry of 
the spectrum due to the parity of the confining potential 
and the interaction potential. Therefore, in order to com- 
pute the NNS distribution we use eigenstates of a given 
parity and spin. This kind of decomposition is a standard 
procedure followed in the analysis of spectral properties 
of quantum systems whenever the Hamiltonian of the 
system possesses a discrete symmetry.Ll After unfolding 
the spectrum the NNS distribution is computed for the 
even parity states. Again, we consider ^ 1000 eigenstates 
of the interacting Hamiltonian (whose eigenenergies are 
lower than the energy of the first transversal mode, for 
compatibility with the Born approximation). 

Since the singlet is the ground state of the two-electron 
system, we first concentrate on the subspace of spatial 
wave functions that are symmetric under particle ex- 
change. The interaction affects very clearly the spec- 
trum, resulting in a strong level repulsion: the NNS 
distribution is in accordance with the predictions of 
RMTJIj As a consequence, the obtained NNS distri- 
bution Pis{s) (histogram shown as a thick soHd hae. 
in Fig. 1) is well described by the Wigner surmiseEHl 
P\y{s) = ^TT s exp(— 7rs^/4), shown for comparison as a 
thin solid line. 

For the triplet states, due to the antisymmetry of the 



spatial wave functions, it is reasonable to expect that 
the tendency of the two electrons to avoid each other re- 
sults in a weaker level mixing. Nevertheless, although 
some differences with the singlet case would appear for 
other statistical measures (that are not possible to per- 
form with the number of levels at hand) , those differences 
are not qualitatively visible on the computed NNS distri- 
bution Pit{s), shown as a dashed line in Fig. 1. Again 
the obtained histogram fits quite well with the predic- 
tions of RMT. 




FIG. 1. NNS distribution obtained for the interacting 
two-electron system described in the text. The solid thick 
line (dashed line) is Pis{s) (P/t(s)). The Wigner surmise 
Pw{s) is plotted as a thin solid line in order to compare it to 
the obtained distributions. Inset: NNS distribution Pni{s) 
for the non-interacting two-electron systems described in the 
text (thick line) togheter with the Poisson distribution Pp{s) 
(thin line). 

We now turn to the dynamics of the classical counter- 
part of the two-electron quantum dot system. We con- 
sider the classical interaction potential given by T4z(|2:i — 
Z2I) = Oij \J{(P -\- \z\ — 2;2p), where the parameters ot and 
d have been obtained from the best fit to the Coulomb 
interaction Ficd^^i ~ Z2\)-, Eq. (||). (The fit is very good 
at all distances, down to the resolution of the spatial 
grid used in our numerical calculations.) The classical 
single particle confining potential Yciz) is a square well 
of length L, and we restrict the analysis to bounded mo- 
tion within this box. In this situation, the effect of the 
confining potential is to reflect elastically the particles 
off the boundaries in each bounce, breaking the transla- 
tional symmetry of the problem. As a consequence, the 
center of mass momentum is not preserved. Nevertheless, 
in the absence of interaction each single particle energy 
is a constant of motion, and therefore the classical prob- 



lem is integrable. On the other hand, the inclusion of the 
Coulomb interaction breaks the conservation of the single 
particle energies and can induce an irregular dymamics 
in the confined system. 

For a total energy E — E* = a/d there is a separatrix 
in the classical dynamics. That is, for E < E* the par- 
ticles never cross each other and the sign of the relative 
coordinate Z2 — zi never changes, while for E > E* it can 
change. 

Denoting e = E/E* we write 



Vc{z[) + Vc{-A) 



d* 



where we have defined 
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(4) 



(5) 



with i = 1,2. In this way, for a given value of the 
reescaled energy e, the classical dynamics depends only 
on the parameter d* . Taking into account that L = 800A 
and the best fit with the quantum Coulomb term Vm 
gives d = 8A, we obtain d* = 0.01. 

In Fig. 2(a) we show for e = 0.9 the Poincare surface of 
section Wj vs. Zj, for the motion of one of the particles, 
taken at times when the other particle bounces off the 
left boundary of the well (the topology of the Poincare 
section does not depend on which particle is selected). 
The motion is chaotic over most of the accesible phase 
space for the given energy shell. Fig. 3(a) shows another 
Poincare section for e = 16. Again, except for small 
regions of regular motion, the dynamic is fully chaotic. 

Although the two Poincare sections look qualitatively 
similar, the trajectories in the plane z'2 vs. z[ are quite 
different as can be seen from Figs. 2(b) and 3(b). Figure 
2(b) (3(b)) shows for e = 0.9 (e = 16) a piece of a tra- 
jectory in the z'2 vs. z[ plane corresponding to an initial 
condition in the chaotic region. In Fig. 2(b) the trajec- 
tory never crosses the straight line defined by z'2 = z'^, 
because for e — 0.9 is always z'2 > z'^. In Ref. O the 
authors perform a classical analysis of the emergence of 
chaos due to the inclusion of an interparticle screened 
Coulomb interaction in an infinite well. The classical 
motion of such a system is qualitatively similar to our 
classical model only for e < 1. 

From the « 1000 eigenstates employed to compute the 
NNS distribution displayed in Fig. 1, only the lowest 
5 % have eigenenergies that correspond to e < 1. The 
eigenenergies of the remaining states correspond to val- 
ues of e ranging from 1 to 16, for which, as we have 
shown, the classical dynamics is chaotic over most part 
of the energy shell. Therefore the NNS distribution com- 
puted from these states results in a remarkable quantum 
signature of the underlying classical chaotic dynamics. 

For energies e >> 1 the classical regular regions 
in the Poincare maps should become predominant over 



the chaotic ones, because in such a limit the Coulomb 
term can be considered as a perturbation of the non- 
interacting two-electron Hamiltonian. Nevertheless, val- 
ues of e >> 1 are not realistic for quantum wells describ- 
ing semiconductors nanoestructures. 
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FIG. 2. (a) Poincare surface of section 112 vs. Z2 for a 
rescaled energy e = 0.9. (b) For the same value of e, a piece of 
a trajectory in the z'2 vs. zj plane corresponding to an initial 
condition in the chaotic region. 
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FIG. 3. (a) Poincare surface of section v'2 vs. z'2 for a 
rescaled energy e = 16. (b) For the same value of e, a piece of 
a trajectory in the z'2 vs. z'l plane corresponding to an initial 
condition in the chaotic region. 

In conclusion, we show that the Coulomb interaction 



is responsible for the chaotic dynamics in a quasi-one- 
dimensional two-electron-quantum dot. This is the first 
study of a realistic few body system where the emer- 
gence of chaos due to interparticle interactions is un- 
mistably demostrated through the analyses of both its 
quantum spectral properties and the dynamics of its clas- 
sical counterpart. We believe that the present results 
also put serious constraints to some models of semicon- 
ductors nanoestructures in which the interaction among 
particles is modeled, for a finite number of particles, as a 
capacitive term in the form of a constant interaction (CI). 
The inclusion of the CI gives an interacting Hamiltonian 
whose spectral properties are those of the non- interacting 
one. In other words, a Poisson NNS distribution will re- 
main Poissonian after considering the interaction in the 
CI model. 
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